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ABSTRACT 

The Rapid Burster (MXB 1730-335) is a unique object, showing both type I and type II X-ray 
bursts. A type I burst of the Rapid Burster was observed with Swift/XRT on 2009 March 5, showing 
photospheric radius expansion for the first time in this source. We report here on the mass and 
radius determination from this photospheric radius expansion burst using a Bayesian approach. After 
marginalization over the likely distance of the system (5.8-10 kpc) we obtain M=l.l±0.3 M and 
R=9.6±1.5 km (la uncertainties) for the compact object, ruling out the stiffest equations of state for 
the neutron star. We study the sensitivity of the results to the distance, the color correction factor, 
and the hydrogen mass fraction in the envelope. We find that only the distance plays a crucial role. 
Subject headings: X-ray, Neutron stars, X-ray binaries, X-ray bursts 



X 

Ok 
6 



> 

m 

O 

(N 



1. INTRODUCTION 



X-ra y bursts were discovered in 1975 bv lGrindlay et all 
(|1976f ) in observations performed with the Astronomi- 
cal Netherlands Satellite of a previously known X-ray 
source, 4U 1820-30. Similar episodes had been observed 
in 196 9 in the Norma co nstellation, with two Vela satel- 
lites (|Belian et alj fl972). but the observed feature was 
related to an accretion event and not identified as a new 
type of source until 1976 (Belian et al. 1976; see also 
Kuulkers et al. 2009 for a reanalysis of that event). The 
first discovered burst (in Cen X-4) is still the nearest- 
known X-ray burst source (at 1 kpc) and the brightest 
X-ray burst ever recorded (-50 Crab0). These pioneer- 
ing discoveries were soon followed by the identification 
of three additional bursters, one of which being the enig- 
matic Rapid Burster (MXB 1730-335, Lewin 1976). To 
date, more tha n 90 Galactic X -ray burst sources have 
b een identified dLiu et alj|2007ft . 

Ivan Horn fc Hansen! (|1974h and lHansen fc van Horn! 
were the first to suggest that nuclear burning on 
the surface of neutron stars could be unstable. Once 
X-ray burst sources were successfully identified in na- 
ture, Maraschi & Cavaliere (1976), and independently, 
Woosley & Taam (1976), proposed that X-ray bursters 
could be powered by thermonuclear runaways on the sur- 
face of accreting neutron stars. However, it was soon re- 
alized that the quick succession of flashes exhibited by 
the Rapid Burster (with recurrent times as short as 10 
sec), did not match the general pattern shown by the 
other bursting sources. A classi fication of type I and type 
II bursts was then established (Hoff man et al.lfl978l ). the 
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former associated with thermonuclear flashes, the later 
linked to accretion instabilities. 

The Rapid Burster is located in the Galactic bulge 
globular cluster Liller 1 (|LewinllT976[ ). at a distance be- 
tween 5.8 and 10 kpc from the Earth (jOrtolani et al.l 
I2007f). A radio transient counterpart was detected 



( Moore et alj |2000) but no IR or optical co unterpart has 
been identified so far ()Homer et al.|[200ll ). It is one of 
the only two known sources showing type II bursts, to- 
gether with the Bursting Pulsar (GRO 1744-28), but the 
Rapid Burster is unique in that it exhibits both type I 
and type II X-ray bursts (for a review see Lewin et al. 
1995). Type I (thermonuclear) bursts display a spectral 
softening during the burst tail, reflecting a decrease in 
the effective temperature and thus the cooling of the neu- 
tron star atmosphere. In contrast, type II bursts exhibit 
a constant temperature for the duration of the burst, 
which can range from —2 s to —700 s, with intervals be- 
tween bursts from —7 s to —1 hr. Over long time-scales, 
as daily monitored by the All Sky Monitor on board the 
Rossi X-ray Time Explorer, the Rapid Burster usually 
radiates at low X-ray flux (< 10 mCrab) indicating a 
low mass accretion rate. Every 100-200 days it enters 
into an outburst Q, reflected in a higher persistent emis- 
sion that indicates a higher mass accretion rate, frequent 
type II X-ray bursts, and some type I X-ray bursts. Out- 
bursts start with a fast rise up to 100 mCrab, followed 
by a decay over two to four weeks. We triggered Swift 
Target of Opportunity observations to follow the burst- 
ing behavior during the January-March 2009 outburst. 
Type II bursts are present in all our observations, but 
their analysis is not the scope of the present paper. 

Here we report on one type I (thermonuclear) burst 
on 2009 March 5 with expansion of the photosphere (as 
shown in section 3), which allows us to constrain the 
mass and radius of the neutron star in the Rapid Burster. 
I Galloway et al.l (|2008f) analyzed all 66 type I X-ray bursts 
of the Rapid Burster available in the public archive of 
RXTE, and found none exhibiting photospheric expan- 

7 Note that in the X-ray binary context, outbursts refer to a 
period of higher flux due to an increased accretion rate, not to be 
confused with explosive, thermonuclear events like nova outbursts. 
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sion. The event analyzed in the present work is thus 
the first one that makes it possible to constrain the mass 
and radius of the neutron star in the Rapid Burster. In 
section 2 we present the observations and the analysis 
methods; section 3 shows the spectral analysis of the pho- 
tospheric radius expansion event (PRE); the determina- 
tion of the mass and radius of the neutron star with their 
uncertainties is presented in section 4; we discuss the re- 
sults and study their sensitivity to the distance, color 
correction factor and hydrogen mass fraction in section 
5; section 6 contains a summary and the conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

Daily Swift X-ray Telescope (XRT, Burrows et al 2005) 
monitoring observations were performed between Febru- 
ary 2 6 and March 5, a ccumulating more than 20 ks of 
data (jSala et al.1 12009) . All observations show intense 
type II bursting activity, with burst intervals between 
25 and 200 seconds. Due to the high average flux of 
the source, all XRT observations were taken in Window 
Timing (WT) mode, which minimizes the pile-up at the 
brightest peaks of the bursts. In WT mode data are free 
of pile- up for count-rates below 100 cts/s. As explained 
below, our results are based on time resolved spectral 
analysis, with the two spectra of interest having count- 
rates of 62 and 25 cts/s, well below the threshold for 
pile-up in WT mode. Swift XRT software was used for 
data reduction (xrtpipeline) and the HEASOFT tool xs- 
elect for light-curve and spectra extraction. For each of 
the spectra, individual exposure maps and ancillary files 
were created, and the XRT response matrix v013 was 
used for all spectral fits. XSPEC 12.6 was used for spec- 
tral analysis. Due to the small number of counts in our 
spectra, we use the C-statistic (|Cashlll979l ) in all fits. 

Most of the detected type II bursts peak at more 
than 100 cts/s in the first Swift/XRT (WT) sequences 
(Feb. 26), which implies fluxes of about 10 -8 erg s^ 1 
cm~ 2 in the range 1-10 keV assuming blackbody emis- 
sion (Lp ea k = (1— 1.6) x 10 38 erg/s at 10 kpc). The peak 
intensity decreases in the last Swift observations, with 
most bursts showing maximum intensities of less than 
60 cts/s in XRT (WT) on 2009 March 5. 

Figure Q] shows the Swift/XRT light-curve of the Rapid 
Burster during the first observation sequence in 2009 
March 5 (OBS ID: 00031360008). A bright type I X- 
ray burst is observed among six other fainter, type II 
bursts. For time resolved spectral analysis of the cooling 
tail of this type I burst, spectra have been extracted for 
a total of 12 time bins, with a variable duration depend- 
ing on the rate, ranging from 0.9 to 20 seconds (Figure 
[2]). The time bins have been carefully selected to follow 
the variations of flux during the burst while keeping the 
maximum possible statistical quality. 

For each time bin, a source spectrum is extracted, se- 
lecting the events in a circular region with a radius of 
20 pixels (~47 arcsec) centered at the source position. 
It is common in burst studies to subtract off the con- 
tinuum, persistent emission from the burst spectrum, 
but this is known to be a source of potential errors 
van Paradijs fc Lewinl 11986( 1 . In our case, the persis- 
tent emission during the orbit containing the type I burst 
(corresponding to all time intervals free of type II bursts) 
is at least two orders of magnitude fainter than the burst 
spectrum, even for the faintest time bins considered for 
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Fig. 1.— Swift/XRT light curve (0.2-10keV) on 2009 March 
5, during the observation sequence containing the PRE burst. 
Six type II burst can be seen in addition to the PRE type I 
burst. The time bin is Is, the start time of the observation is 
14895 at 13:27:13.216 in truncated Julian Days (TJD, Julian Day 
- 2440000.5). 

the mass and radius determination in Section 4. To check 
the possible effect of the persistent emission in our re- 
sults, we performed the spectral analysis for the faintest 
spectra of the type I burst with and without persistent 
emission subtraction. The resulting spectral fit param- 
eters are identical within less than 1%, much smaller 
than our final uncertainties, confirming that the persis- 
tent emission does not affect our results. The systematic 
uncertainty in the Swift/XRT absolut e calibration, con- 
sidered to be 10% (jGodet et al.ll2009T ). is also negligible 
compared to our final uncertainties and does not affect 
our results. 

3. SPECTRAL ANALYSIS 

All burst spectra were fit in the energy range 1-10 
keV with a blackbody model, modified with interstellar 
absorption by the Tubingen-Boulder absorption model 
([Wilms et al.l l2000). This model is used to calculate the 
cross section for X-ray absorption as the sum of the cross 
sections due to the gas-phase, the grain-phase, and the 
molecules in the interstellar medium. The hydrogen den- 
sity column is kept co nstant to its average inte rstellar 
value, 1.7 x 10 22 cm~ 2 (lDickev fc Lockmanl 119901 ) Not e 
that using the relation of iPredehl fc Schmittl (|1995f ). 
this column density corresponds to an extinction of 
E(B-V)=3, as obser ve d for the Liller 1 glob ular cluster 
(|Barbuv et al.l I1998Q . iMarshall et~aT1 (|200lD found the 
same absorbing hydrogen column density to the Rapid 
Burster with Chandra/HETG spectra of type II X-ray 
bursts, indicating that there is no source intrinsic absorp- 
tion. The blackbody model used for the fit of the cooling 
tail spectrum is parametrized with the effective temper- 
ature and a normalization constant proportional to the 
emitting area, more precisely = {R^^ / D wkpc ) 2 

(with -R^ dy being the radius in km and Diofcpc the dis- 
tance to the source in units of 10 kpc). The evolution 
of the bolometric flux, the normalization parameter Aqo, 
and the blackbody temperature obtained from the spec- 
tral fits is shown in Figure [2] 

Three clear phases are visible in the type I burst light 
curve and spectral parameters evolution, clearly indicat- 
ing that photospheric radius expansion (PRE) is occur- 
ring. First, the expansion of the photosphere occurs dur- 
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Fig. 2. — Time resolved spectral analysis: evolution of the bolo- 
metric flux and the spectral parameters during the PRE burst. 

ing the initial rise in flux, lasting ~ 3 seconds. Second, 
after maximum expansion, the contraction of the pho- 
tosphere follows, lasting ~ 7 seconds. The contraction 
finishes at the so called touchdown of the photosphere 
(occurring at bins 6-7 in Figure[2J 8-12 seconds after the 
start of the burst). At this moment, the photospheric 
radius returns to the neutron star surface and the black- 
body radius reaches its minimum value and the color 
temperature its maximum. Finally, the cooling down of 
the envelope is observed (time bins 8-11, 13-54 seconds 
after the start of the burst). The twelfth and last bin 
analyzed and shown in Figure [D does clearly show that 
neither the flux nor the temperature are decreasing any- 
more, so we consider the cooling tail to finish 54 seconds 
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Fig. 3. — Confidence contours at 68%, 90% and 99% confidence 
level for the spectral parameters at touchdown. Normalization is 
in units of L39/D 2 0kpc , with L39 being the bolometric luminosity 
in units of 10 39 erg/s and Di 0kpc the distance in units of 10 kpc. 
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Fig. 4. — Confidence contours at 68%, 90% and 99% confidence 
level for the spectral fit parameters during the cooling tail. Nor- 
malization is in units of {f^°f v / D wkvc ) 2 , where R*£^ y is the 
apparent radius in km and D\okpc the distance to the source in 
units of 10 kpc. Numbers by the contours indicate the time bin 
(as shown in Fig [2j . Diagonal lines corresponding to constant lu- 
minosity are shown for luminosities 10 37 and 10 3s erg/s. 

after the start of the burst. 

We devote some attention to quantify the radius ex- 
pansion. A straightforward calculation of the x 2 value 
with respect to the average radius would not be cor- 
rect, because due to the low statistics of the data, the 
probability density function (PDF) of the spectral pa- 
rameters in each of the time bins follows a Poisson dis- 
tribution and differs clearly from a normal distribution. 
In other words, error bars are not symmetric, as shown 
for the touchdown and cooling tail bins in Figures [5] and 
[H To deal properly with the PDF of each of the bins 
of the radius evolution and evaluate the radius expan- 
sion, we have taken into account the calculated PDF as- 
sociated to each bin and fit each of them by a Poisson 
distribution. We have then performed a model selec- 
tio n using the Bayesia n algorithm FABADAEI proposed 
in iPardo et al.l (|2011al) . We first fit the radius evolu- 

http:/ /gem. upc.edu/members/luis-carlos/luis-carlos- 
pardo/bayesiano 
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tion shown in Figure [5] with a constant radius (more pre- 
cisely, with ( R^ dy / D-\uk-nr) 2 = 70), redefining \ 2 as m 
IPardo et al.l (|20llbD to take into account the asymmetry 
of the PDF describing the data. The x 2 obtained mod- 
eling the data by a constant radius is 16.5 (for 5 degrees 
of freedom), with a chance probability of only 0.55%. In 
contrast, a model consisting of an expansion and con- 
traction of the radius provides x 2 =3.5 (for 3 degrees of 
freedom). Therefore, the expansion is clearely favored in 
comparison to a constant radius. 

The bolometric luminosity is almost constant during 
the envelope expansion (0-3 s) and contraction (3-12 s), 
and decreases during the cooling tail (11-53 s). Expan- 
sion of the photosphere is caused by the luminosity reach- 
ing the Eddington limit. The neutron star mass can then 
be determined from the observed bolometric luminosity, 
assuming that the touchdown (i.e. maximum) luminosity 
corresponds to the Eddington limit as seen by a distant 
observer. After the touchdown, the emitting area is as- 
sumed to be the whole neutron star surface, and from 
the spectral fits during the cooling tail, the neutron star 
radius can be inferred. For convenience in the analysis, 
since we are interested in the cooling emitting area and 
the touchdown flux, the blackbody model used for the 
cooling tail is parametrized with the normalization pro- 
portional to the emitting area {Aoo = (i?^ dy /Z)iofepc) 2 ) 
as explained above for all time bins), while for the touch- 
down spectrum we use a normalization proportional to 
the bolometric flux as L39 / D\ Qkpc , where L 3g is the bolo- 
metric luminosity in units of 10 39 erg/s. 

For the determination of the touchdown flux, we merge 
time-bins 6 and 7 of Figure[5]to obtain a higher signal-to- 
noise spectrum. Figures [3] and @] show the confidence con- 
tours of the spectral parameters obtained for the touch- 
down spectrum and the spectra of the four time bins 
of the cooling tail. Since the count-rate is decreasing 
during the cooling tail, the parameters become less well 
constrained for the later time-bins (labeled as 9, 10 and 
11 in Figure 0J. For the determination of the cooling 
emitting area, we consider the range of values for the 

Aoo = (R b km dy I Dwkpc) 2 obtained for the first bin of the 
cooling tail (bin 8), which has the brightest spectrum and 
thus the better constrained emitting area. Note that all 
later bins, while having larger uncertainty ranges, are 
compatible with the results obtained from bin 8. 

4. MASS AND RADIUS OF THE NEUTRON STAR IN THE 
RAPID BURSTER 

The time resolved spectral analysis of X-ray data of 
a PRE burst initially provides only the best-fit black- 
body parameters. But there are a number of factors that 
play a role in the determination of the actual neutron 
star radius and mass. Firstly, the neutron star spec- 
trum is known to differ from a perfect blackbody, with 
the bursting neutron star atmosphere emitting a harder 
spectrum than a blackbody due to electron scattering 
in the hot neutron star atmosphere and th e resulting 
suppression of emissivity ([van Paradiisl 11982) . The ac- 
tual effective temperature can be related to the color or 
blackbody temperature obtained from the spectral fit by 
the color correction factor, f c = Tbbody/T e ff. A num- 
ber of recent works have been devoted to the modeling 
of the neutron star atmosphere spectra and the deter- 
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Fig. 5. — Probability density function of the touchdown flux. 
Normalization is in units of L39 / D^ okpc , where L39 is the bolo- 
metric luminosity in units of 10 39 erg/s. A Gaussian is shown in 
red, dashed line for comparison. 
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Fig. 6. — Probability density function of the cooling area. Nor- 
malization units are (fl^^/Diofcpc) 2 - A Gaussian is shown in 
red, dashed line for comparison. 

mination of t (IMadei et al.ll2004t IMaiczvna et al.H2005t 
ISuleimanov et aLl l2011). The color correction factor f c 
will affect both the temperature and the emitting area, 
since this is obtained from the observed flux and black- 
body temperature. A second factor affecting the radius 
and mass determination of the neutron star is the dis- 
tance D, needed for the luminosity determination. The 
hydrogen mass fraction X plays its role in the Eddington 
limit for the luminosity through the electron scattering 
opacity, K e — 0.2(1 + X)cm 2 g~ 1 . Finally, the high sur- 
face gravity of the neutron star causes light deflection, 
affecting the observed flux and emitting area. In a PRE 
burst, the gravitational rcdshift changes with the expan- 
sion of the photosphere and special caution is needed in 
the selection of the moment for the determination of the 
bolometric flux. The best moment for the determination 
of the bolometric flux is at touchdown, when the luminos- 
ity is still at the Eddington limit, but the photospheric 
radius is already at the neutron star surface, with the 
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redshift then determined by the neutron star radius. 

Combining the above considerations, the observed 
touchdown flux F^ and the apparent angular emitting 
area Aoo depend on the actual neutron star radius R and 
mass M as 

GMc f 2GM\ 1/2 
00 ~ 0.2(1 + X)D? \ (1) 

and 

Aoo ~ D*f* \ L Re* ) [2) 

With the uncertainty ranges of our observed touch- 
down flux and cooling emitting area as shown in the 
confidence contour plots in Figures [3] and |H an accu- 
rate treatment of the error ranges is needed. Here we 
follow the Bayesian approach introduced in lOzell (|2006T > 
and lOzel et al.l ()2009f ) in the first two of a series of papers 
by these authors on the determination of the mass and 
radii for several bursting neutron stars. A probability 
density function (PDF) is assigned to each of the observ- 
able quantities (touchdown flux, cooling area, distance), 
designated by P(F 00 )dF 00 , P(A 00 )dA 00 , P(D)dD, as well 
as the color correction factor and the hydrogen mass frac- 
tion, P(f c )df c and P(X)dX. Since these are all indepen- 
dently determined parameters, the combined probability 
for given values of Foo, A^, D, f c and X is directly the 
product of the individual probabilities. This combined 
PDF is then converted to the PDF of the neutron star 
mass and radius using the change of variables as given in 
equations (1) and (2). 

We assign uniform probabilities to the distance, be- 
tween 5.8 and 10 kpc (Ortolani et al. 2007, sec how- 
ever next section); to the mass hydrogen fraction, be- 
tween and 0.7 (solar); and to the color correction 
factor, between 1.3 and 1.4. This last choice is based 
on the values obtained from theor etical modeling of 
bursting neutron star atmospheres (|Madei et all 120041 : 
iMajczyna et al.l 120051 : ISuleimanov et al.l 120 111 ) and the 
fact that the ratio L/L^c/ at the time we measure the 
emitting area is in the range 0.03-0.4. 

For the observable parameters F^ and A^, instead 
of assigning a Gaussian distribution with the values of 
the spectral fit, we determine the PDF of the observed 
touchdown flux, P(F 00 )dF 00 , and of the cooling area, 
P(A 00 )dA 00 , directly from the blackbody model fit to the 
Swift/XRT spectra by running a Monte Carlo Markov 
chain (MCMC) to each of the spectral fits (touchdown 
and cooling). This choice is based on the fact that the 
results of the MCMC show that both P(A 00 )dA 00 and 
P(F 00 )dF 00 are asymmetric functions and differ consid- 
erably from a Gaussian distribution, as shown in Figures 
Eland© 

To compute the final PDF on the M-R plane, we pro- 
ceed as follows: for each M-R pair we calculate the cor- 
responding values of A^ and F^ for given values of 
D, f c and X within the range considered, with equa- 
tions (1) and (2); we then determine the P(A 00 )dA 00 
and P(F 00 )dF 00 corresponding to these values from our 
observational results and calculate the combined PDF. 
We marginalize over (integrate over all values of) D, f c 
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Fig. 7. — Mass and radius of the neutron star in the Rapid 
Burster as determined from the type I X-ray burst with Photo- 
spheric Radius Expansion observed by Swift/XRT on 2009 March 
5. The color map shows the probability density function after 
marginalization over the whole range of the distance (5.8-10 kpc), 
the color correction factor (1.3—1.4) and the hydrogen mass fraction 
(0—0.7). Confidence contours at la and 2a are shown for the dis- 
tance ranges 5.8—6.5 kpc, 7.5-8.3 kpc, and 8—10 kpc, corresponding 
to different distance determinations in Ortolani et al. (2007) and 
Valcnti et al. (2010). Mass-radius relations for several equations 
of state are shown as well, with relations labeled with MS0, AP3, 
UU, FPS and A corresponding to equations of state of neutron 
stars without condensates, M, GS1 and GS2 to equations of state 
with condensates, and AF0, SSI and SS2 to strange quark matter 
stars (see Lattimer and Prakash 2001 for more details). The upper 
left part of the plane (left of the solid line) is excluded by causality. 

and X to obtain the probability density function on the 
M-R plane as shown by the color map in Figure [7l In the 
next section we test the robustness of the results to vari- 
ations of the distance, color correction factor and mass 
hydrogen fraction. 

Finally, we obtain the PDF of the neutron star mass by 
marginalizing over the radius, and the PDF of the radius 
by marginalizing over the mass. The final radius PDF is 
well fit with a lognormal distribution with median 9.6 km 
and a variance of 2.5 km (a Gaussian fit provides a similar 
mean 9.6 km and <7=1.5 km). The final PDF of the mass 
is acceptably represented by a Gaussian distribution with 
mean 1.1 M and er=0.3 M©. 

5. DISCUSSION 

Time resolved spectroscopy of the type I X-ray burst 
from the Rapid Burster with photospheric radius expan- 
sion allows us to constrain the mass and radius of the 
accreting neutron star by determining the probability 
density function on the mass-radius plane, as shown in 
Figure [7] For the accreting neutron star in the Rapid 
Burster we obtain M=l.l±0.3 M and R=9.6±1.5 km 
(la uncertainties) for the assumed distance range 5.8- 
10 kpc. This indicates that the accreting object in the 
Rapid Burster is compatible with both normal matter 
and strange quark matter equations of state, but the 
stiffest equations of state can be ruled out. 

This is, to our knowledge, the first detected type I 
X-ray burst of th e Rapid Burster show ing photospheric 
radius expansion. iGallowav et al.l (|2008f) analyzed all the 
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Fig. 8. — Sensitivity of the final mass and radius of the neutron 
star to the color correction factor f c (upper panel) and the hydrogen 
mass fraction X (lower panel). Confidence contours at la and 2<r 
are shown for f c = 1 . 3 (solid lines) and f c =1.4 (dotted lines), and 
for X=0 (solid lines) and X=l (dotted lines), while the color scale 
shows the probability distribution function as in Figure [7] for f c 
marginalized over 1.3—1.4 and X over 0—0.7. 

observations of the Rapid Burster available in the RXTE 
archive and identified 66 type I X-ray bursts. While 
some of those bursts show a flat topped light curve, none 
of them follows a temperature and radius evolution like 
the one shown in our Figure [2] In addition, the peak 
flux of the X-ray burst analyzed in the present work is 
about a factor 3 higher than the peak flux of the bright- 
est type I X-ra y burst s of the Rapid Burster shown by 
I Galloway et al.l (|2008l) . supporting that in the present 
case the Eddington limit of the luminosity was reached. 

We investigate now in detail the effect of these choices 
and the sensitivity of the final result to these parame- 
ters. We con sider first the e f fect o f the distance on the 
final results. lOrtolani et al.l (|2007fl determined the dis- 
tance to Liller 1 based on HST NICMOS F110W and 
F160W band photometry and by using the location of the 
horizontal branch in the color-magnitude diagram com- 
pared to other globular clusters with known distances. 
The F110W and F160W filters are different than the 
ground-based equivalent J and H band filters in that 
they are wider and extend further to the blue than the 
ground-based filters. A careful transformation of NIC- 
MOS photometry into a ground-based system is there- 
fore necessary. This is particularly true for observations 
of Liller 1, whic h suffer from some fo reground extinction 
of E(B-V)= 3.0 dBarbuv et al.lll998l) . Unfortunately, as 
indicated in lOrtolani et al.l (|2007fl . the choice of the color 
transformation is of importance, in particular for the J-H 
c olor term. 

lOrtolani et al.l (|2007| ) determined the distance to Liller 
1 assuming metallicities [Fe/H]=-0.1 and [Fe/H]=-0.6, 
and for three different photometric calibrations, report- 
ing a total of six dista nce estimates in the range 5.8- 
10.0 kpc. More recently. I Valenti et al.l(|2010D re-observed 



Liller 1 in the near-IR using a standard ground-based 
filter set where a transformation to the 2MASS system 
is well understood. They found a distance modulus of 
14.48 mag (d=7.87 kpc) with an uncertainty of 0.2 mag 
(0.4 kpc), in the center of the ran ge of values for the 
distance from lOrtolani et "all (|2007t ). 

The results for the neutron star mass and radius deter- 
mined for the distance ranges 5.8- 6.5 kpc (Stephens e t 
al. 2000 calibration), 7.9(±0.4) kpc (| Valenti et al.ll20"l0l) . 
and 8-10 kpc (other calibrations) are shown in Figure [2 
The mass and radius o btained for th e most recent dis- 
tance determination by iValenti et al.l (|2010T ) are consis- 
tent with the values fou nd taking into a ccoun t the whole 
uncertainty range from lOrtolani et "all ()2007l ). 

The hydrogen mass fraction X plays a role in the Ed- 
dington limit for the luminosity (equation 1). To ob- 
tain the results shown in Figure [7] we have integrated the 
PDF for all possible values of X between zero and so- 
lar (X=0.7) assuming solar matter is accreted, and then 
metallicity subsequently increased due to the nucleosyn- 
thesis during the thermonuclear X-ray burst. To test the 
sensitivity of the results on the mass-radius plane, we 
have calculated the probability distribution function for 
two fixed, extreme values of X, X=0 (no hydrogen) and 
X=l (pure hydrogen, even if this is an unrealistic case). 
Results are shown in Figure [5] A high hydrogen mass 
fraction corresponds to higher masses and radii. Never- 
theless, the fact that the hydrogen fraction plays its role 
in the determination of the Eddington limit is reflected 
in a higher sensitivity of the final mass to the value of X. 
Since all possible values for the hydrogen mass fraction 
predict physically plausible results for the mass and ra- 
dius, this test provides no constraints on X for the Rapid 
Burster. 

Finally, the result of a similar test on the color cor- 
rection factor is also shown in Figure [8] In this case, 
a higher color correction factor corresponds to higher 
mass and radius. Since the color correction factor plays 
a role in the determination of the neutron star radius 
from the emitting area (equation 2), the radius is the 
parameter mor e affected. Note that works on atmo- 
sphere models dMadei et al.ll200l iMaiczvna et al.ll2005l : 
iSuleimanov et al.l 1201 ID calculate the color correction 
factors only for metallicities of solar values or smaller. 
The smallest values for the color correction factor are 
obtained for solar abundances, i.e, the highest metallic- 
ity considered in those works. The small color correction 
factor is due to the absorption edges of Fe, which move 
the neutron star hardened spectrum back towards lower 
energies, closer to the blackbody. However, in Edding- 
ton limited bursts, the photosphere of the neutron star 
is expected to show ashes from H/He burning, produced 
deep in the envelope (Weinberg, Bildsten & Schatz 2006, 
in 't Zand & Weinberg 2010). This is in agreement with 
results from hydrodynamic simulations, which show that 
convective mixing reaches the surface only when the flux 
reaches the Eddington limit (Woosley et al. 2006, Fiskcr 
et al. 2008, Jose et al. 2010). This would imply more 
absorption features in the final emerging spectrum and 
a smaller color correction factor, thus our resulting PDF 
would move down to smaller masses and radii. So we can 
consider that the color correction factors taken into ac- 
count in the present work provide, in the worst case, an 
upper limit for the mass and radius of the Rapid Burster. 
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6. SUMMARY 

The Rapid Burster showed a type I X-ray burst with 
photosphcric radius expansion on 2009 March 5. The 
analysis of Swift/XRT time- resolved spectra of that burst 
allows us to constrain the mass and radius of the neu- 
tron star in the Rapid Burster. After marginalization 
over the whole uncertainty range for the distance (5.8- 
10 kpc), we obtain M=l.l±0.3 M and R=9.6±1.5 km 
(lcr uncertainties), compatible with both normal mat- 
ter and strange quark matter equations of state for the 
accreting object in the Rapid Burster. We check the sen- 
sitivity of the results to the distance, the color correction 
factor, and the hydrogen mass fraction, and we find that 
only the distance plays a role. Most recent distance de- 
terminations locate the Globular Cluster Liller 1 hosting 
the Rapid Burster at the center of the whole uncertainty 
range, confirming the results for the mass and radius 
found above. 

As far as we know, the type I burst shown in the 
present work is the sole example from the Rapid Burster 
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showing photospheric radius expansion. Therefore, fu- 
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new type I bursts with photospheric expansion would 
provide a confirmation of the results presented. We 
have checked that neither the uncertainties in the hy- 
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can change our conclusions. Only the distance is the key 
factor in this case, and thus a better determination of the 
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requires a definitive measurement of the distance to Liller 
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